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We propose the Wigner separability entropy as a measure of complexity of a quantum state. This 
quantity measures the number of terms that effectively contribute to the Schmidt decomposition of 
the Wigner function with respect to a chosen phase space decomposition. We prove that the Wigner 
separability entropy is equal to the operator space entanglement entropy, measuring entanglement 
in the space of operators, and, for pure states, to twice the entropy of entanglement. The quantum 
to classical correspondence between the Wigner separability entropy and the separability entropy of 
the classical phase space Liouville density is illustrated by means of numerical simulations of chaotic 
maps. In this way, the separability entropy emerges as an extremely broad complexity quantifier in 
both the classical and quantum realms. 
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I. INTRODUCTION 



Measuring complexity in a simple and unified way has 
been a major and long quest in both quantum and clas- 
sical dynamics. While there exists a direct connection 
between chaos and algorithmic complexity of trajectories 
in classical physics [l[ , the problem is particularly elusive 
for quantum mechanics , where the notion of trajectory 
is forbidden by the Heisenberg uncertainty principle and 
complexity can be attributed not only to the lack of inte- 
grability but also to the tensor-product structure of the 
Hilbert space, that is, to entanglement. 

The phase space representation of quantum mechan- 
ics is a very convenient framework to investigate quan- 
tum complexity, in that one can compare classical and 
quantum dynamical evolutions of distributions in phase 
space. In this context, the number of Fourier harmonics 
of the Wigner function has already been used to char- 
acterize the complexity of a quantum state, both for 
single-particle Q and many-body Q quantum dynam- 
ics, in particular to detect quantum phase transitions. 
This complexity measure can be equally applied to clas- 
sical and quantum mechanics, with the Liouville den- 
sity used instead of the Wigner function in the classi- 
cal case. However, this quantity has the disadvantage 
of being basis-dependent. Moreover, knowledge of the 
whole Fourier harmonics spectrum of the Wigner func- 
tion seems in general not necessary to compute expec- 
tation values of physically relevant observables. For in- 
stance, even though the number of harmonics grows ex- 
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ponentially in time for both integrable and non-integrable 
quantum chaotic Ising chains [J|, the resources required 
to simulate both local and extensive observables grow ex- 
ponentially in the chaotic case but only polynomially for 
the integrable model |5|. 

Recently, a new complexity indicator has been intro- 
duced for classical dynamics, the separability entropy 0, 
measuring the logarithm of the effective number of terms 
in the Schmidt (or singular value) decomposition of the 
Liouville density, with respect to an arbitrary phase space 
decomposition. This quantity estimates the minimal 
amount of computational resources required to simulate 
the classical Liouvillian evolution and grows linearly in 
time for dynamics that cannot be efficiently simulated. 
In this paper, we extend this notion of complexity to 
the quantum realm, by defining the Wigner separabil- 
ity entropy as the number of terms that effectively con- 
tribute to the Schmidt decomposition of the Wigner func- 
tion. We prove that such quantity is equal to the op- 
erator space entanglement entropy Q, constructed from 
the Schmidt decomposition of the density operator in the 
space of Hilbert-Schmidt operators and quantifying the 
complexity of time-dependent density-matrix renormal- 
ization group simulations. Furthermore, for pure states 
the Wigner separability entropy is twice the entangle- 
ment entropy. We illustrate the quantum to classical 
correspondence for the separability entropy and its link 
to the entanglement entropy by means of numerical sim- 
ulations of quantum chaotic maps. 

The paper is organized as follows. In Sec. [IT] we define 
the Wigner separability entropy h[W] and prove its rela- 
tion to the operator space entanglement entropy and, for 
pure states, to the entanglement entropy. The dynami- 
cal evolution of h[W] is studied for quantum chaotic map 
models in Sec. Mil We finish with concluding remarks in 
Sec. EH 
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II. WIGNER SEPARABILITY ENTROPY 

Given a system described in a 2<i-dimensional compact 
phase space f2 by the Wigner function W(z) (with the 
normalization constraint J dzW(z) = 1) and an arbi- 
trary phase space decomposition, f2 = Hi © f2 2 , into two 
set of coordinates, z = (x, y), we can write the Schmidt 
(singular value) decomposition of the Wigner function: 



(x)w n (y) 



(1) 



with n € N, {v n } and {w n } orthonormal bases for 
L 2 (f2i) and L 2 (f2 2 ), respectively, and the Schmidt co- 
efficients (singular values) pi > p2 > . . . > satisfying 
J2 n ^n = I dzW 2 (z). We then define the Wigner sepa- 
rability entropy as 



where 



Pn 



Pn 



JdzW 2 (z) 



(2) 



(3) 



The coefficients {p n } satisfy J2 n An = 1 an d are the 
Schmidt coefficients of the singular value decomposition 



oi W = W/ y J dzW 2 (z). That is, W is normalized in 
L 2 (n): JdzW 2 (z) = 1. Note that for pure states 



= (27rn)- d / 2 P^(2). 



(4) 



The Wigner separability entropy h[W] quantifies the log- 
arithm of the number of terms that effectively contribute 
to decomposition (JTJ) and therefore provides a measure of 
separability of the Wigner function with respect to the 
chosen phase space decomposition. 

The main advantage of defining the separability en- 
tropy in phase space by means of the Wigner function is 
that such quantity can be directly translated to classical 
mechanics. The classical analog of the Wigner separa- 
bility entropy is the s-entropy h[p c ] defined in Ref. @, 
where the classical phase space distribution p c {z) is used 
instead of the Wigner function W(z). 

It is interesting to establish a connection between the 
Wigner separability entropy and the operator space en- 
tanglement entropy 7], constructed from the Schmidt 
decomposition of the density operator p acting on the 
Hilbert space H. Since Tr(/5 2 ) < 1, the density operator 
is a Hilbert-Schmidt operator, that is, an operator with 
finite Hilbert-Schmidt norm ||/3||hs = \/Tr(/S 2 ). There- 
fore, p can be interpreted as a vector belonging to the 
Hilbert space B(H) of Hilbert-Schmidt operators, with 
the Hilbert-Schmidt inner product (A,B)ns = Tr(A^B). 
Therefore, given % = Hi <X> H2, the density operator has 
a Schmidt decomposition, 



(5) 



where {cr n } and {f n } are orthonormal (Trfcrj^oy,) = <5 m . n , 
Tr(r^jf n ) = S m<n ) bases for B(Hi) and £?("H 2 ), respec- 
tively, and the Schmidt coefficients p,\ > p,2 > • • • > 
satisfying J2 n Pn = TrQif) = IIpIIhs- The operator space 
entanglement entropy [7j is then given by 



Ihi 



Pn 



\\p\\ 



HS 



(6) 



In what follows, we prove that h[p] = h[W], that is, the 
operator space entanglement entropy equals the Wigner 
separability entropy. This result follows from the fact 
that the Weyl correspondence establishes an isomorphism 
between Hilbert-Schmidt operators and L 2 (fl) functions 
on classical phase space [|[ HJ. Since the density oper- 
ator p is the integral kernel of a unique linear one-to- 
one transformation mapping the operators a n f-> f„ , Vn, 
i.e. Tri((<rJ)i/3) = p n T n , and vice versa for the inverse 
transformation, 1^2 ((7^)2/5) = Pn^m and since the Weyl 
transform of the density operator is the Wigner func- 
tion, it follows that p and W have the same Schmidt 
coefficients, and therefore 



h[p] = h[W\. 



(7) 



In order to explicitly illustrate the above result, we 
consider a density operator p(a\, 0,4, a\, a' d ) written 
in terms of a set of bosonic creation-annihilation opera- 
tors, [dj, aj] — [a J, a] 
Wigner function as 



0, [aj,aji 



Sij, and define the 



W(a, a*) = 



1 



dr/dt]* exp 



77 • a r] a 



j — y ^ ^ 

xTr{pD( V )}, (8) 

where r] — (771, n d ) £ C d and a = (ai, ay) € C d 
are (^-dimensional complex variables, the integration runs 
over the complex ^-planes for i = l,...,d, D is the dis- 
placement operator 



D{-q)=D (771, ...,r)d) = exp 
and j a) are the coherent states 



a 

,i=i 



, (9) 



|a) = \aia 2 ...a d ) = D I -j= j |0), (10) 

with \cti) being an eigenstate of the annihilation operator 
hi, i.e., ai\cti) — ^t|aj), and |0) being the vacuum state 
a,j\0) = 0. Using the singular value decomposition 

P = slttn<7n(a>i, —, ad/2, «j, — ) a d /i) 

n 

®T n {a, d/2+ i, a d , a d/2+v a d ), (11) 

we obtain 

Ti[pD(r])} = ^2p n & n {i ll ,...,i ld/ 2,r]l,...,n* d/2 ) 

n 

Xf n (Vd/2+l, ■■; Vd, Vd/2+V Vd)' ( 12 ) 
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where 

5-« = Th,...,d/2 [&nD(r)i,...,ri d /2)], (13) 
?V» = Tr d/ 2 + i i ... id [f n -D(% /2+1 ,. ..,%)]• (14) 

Finally, we derive 

(a 1 ,...,a d/2 ,a* 1 ,...,a* d/2 ) 

n 

Xw n (a d / 2+ x, a d , a* d/2+1 , a d ),(l5) 



where 



1 



n d h d/2 

' d/2 



x exp 



E 

i=l 



d^dr/* ■ ■ • di] d/2 dr]* /2 



(16) 



/ d^/a+id^ 



x exp 



E 

i=d/2+l 



2+1 



dVddi] d 



(17) 



Eq. (fT5|) is precisely the singular value decomposition 
of the Wigner function (the orthogonality of the Schmidt 
vectors (|16I17I) can be checked straightforwardly) and has 
the same Schmidt coefficients {p n } as the the singular 
value decomposition (1111) of the density operator, so we 
conclude that h[W] = h[p] ©. 

When the density operator p describes a pure state, 
P — IV 1 ) (V'li there exists a simple relation between the 
Wigner separability entropy and the entanglement con- 
tent of the state \if>) G U = 'Hi ® % 2 . The Schmidt 
decomposition of \ip) reads [Tol. Hlj 



E A J 



(18) 



with {<f>j} orthonormal basis for Hi (H2), Ai > 

A2 > . . . > 0, and J^j = 1- On the other hand, we 
can also write the Schmidt decomposition of the operator 

p = ^A 3 A fe |^)(0 fc |®|^)(a|. (19) 

The comparison between (fT9|) and ([5]) implies {/injneN = 
{AjAfe}j,fe e N- Therefore, 

w = -E^« ln ^ = -E A ? A ' ln ( A ' A ') 
= - 2 E A > A ' 



-2S(p 1 ) = -25(p 2 ),(20) 



where pi = Tr 2 (/5) and /5 2 = Tri (/5) are the reduced den- 
sity operators for subsystems 1 and 2 and S 1 is the von 



Neumann entropy. Since for a pure state \ip) von Neu- 
mann entropy of the reduced density matrix quantifies 
the entanglement E of \ip) (ToL ITlj. 



£(|^» = S(pi) = S(p 2 ), 



(21) 



we can conclude that the Wigner separability entropy is 
twice the entanglement entropy E(\ip)): 



h[W] = 2E(\ip)). 



(22) 



For pure states S(p) — 0, and therefore the Wigner 
separability entropy is equal to the quantum mutual in- 
formation 



I(l:2) = S(pi) + S(p 2 )-S(p). 



(23) 



Quantum mutual information measures total correla- 
tions, both of classical and quantum nature, between sub- 
systems 1 and 2 and for pure states classical correlations 
CQip}) are equal to quantum correlations, measured by 

mm [il- 



iii. WIGNER SEPARABILITY ENTROPY FOR 
CHAOTIC MAPS 

In order to illustrate the quantum to classical corre- 
spondence for the separability entropy, we study the evo- 
lution in time of h[W] and of its classical counterpart 
h[p c ] for three quantum maps: (i) the quantum baker's 
map, (ii) the perturbed quantum cat map, and (iii) a 
model consisting of two coupled perturbed cat maps. The 
first model (i) exhibits an atypical failure of the quan- 
tum to classical correspondence due to discontinuity in 
the mapping producing drastic quantum diffraction ef- 
fects. The second model (ii) is typical in that the Wigner 
separability entropy h[W] follows the classical separabil- 
ity entropy h[p c ] up to the Ehrenfest time scale t E , un- 
til Wigner function evolution (whose time derivative is 
given by the Moyal bracket) is well approximated by the 
evolution of the classical Liouville density (whose time 
derivative is given by the Poisson bracket, the first order 
in h expansion of the Moyal bracket), i.e. 

W(t = 0) = p c (t = 0) => h[W(t)} = h[p c (t)} for t < t E . 

(24) 

Note that for chaotic dynamics with Lyapunov exponent 
A, the Ehrenfest time scales as t% ~ (— logfi)/A [13|. The 
third model (iii) illustrates for pure states the connection 
between the separability entropy and the entanglement 
entropy. 

We study the quantum versions of classical chaotic 
maps on the 2-torus [0, 1) x [0, 1). Quantization imposes 
the Planck constant to coincide with the inverse of the 
dimension N of the Hilbert space, i.e., H = 1/{2ttN). The 
baker map has been studied in detail in many works (see 
for example [3, [l5[ ) . The classical baker map is defined 
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by 

(9t+1 ' Pt+l) -\(2g t -l,ip t + i), if ±<g t <l, 

(25) 

where the discrete time t E Z measures the number of 
map iterations. The quantum analogue of the classi- 
cal baker map is a unitary operator acting on an N- 
dimensional Hilbcrt space (assuming N to be even). In 
the positions-representation its matrix reads: 



where q l and p l (i — 1,2) are again taken mod (1), 
and where we have used the same perturbation e(q) 
as before (for both maps). Here, e'(qj,q 2 ) — 
— [K c /(2ir)} sin (2ivql + 2nq 2 ) is the coupling term. 

The quantized version of the 2D perturbed cat map 
is obtained by multiplying a separable (tensor) product 
of quantized ID perturbed cat maps with a simple ex- 
ponentiated coupling matrix (which is diagonal in the q- 
representation) . More explicitly, iV 2 x iV 2 unitary matrix 
M 2D is given by 



-l ( G N /2 



B = G-J ^ N J 



N \ G 



N/2 



(26) 



where Gn is the TV-dimensional antiperiodic Fourier ma- 
trix: 



1 



(G N )ij = -== e -»"WW+W/» , (27) 
v N 



(28) 



with < l,j < N- 1. 

The classical perturbed cat map reads fl( 



Qt+i 
Pt+i 



= M 



1i 
Pt 



where q and p are taken mod (1), and where e(g*) = 
— [K/(2ir)\ sin (2irq t ). Throughout this paper we have 
used 



M 



2 1 

3 2 



(29) 



satisfying the symplectic condition detAi = 1, and 
K = 0.5. The perturbed quantum cat map in the q- 
representation is given by the N x N unitary matrix M 
whose elements read 



Mij = A exp 
where 



NMv. 



-(Mnl 2 -2lj+M22j 2 )+F 



(30) 



A = [1/(\NM 12 )} 1/2 , 

F = [iK N/ '(27T)] cos (2-kI/N). 

Finally, we have considered a 2D model (in 4D phase 
space (q 1 , q 2 , p 1 , p 2 ) € [0, l) 4 ) consisting of two coupled 
perturbed cat maps. In order to do this in a symplectic 
way we have used a symmetric coupling in the position 
coordinates [13] ■ In this case the classical map is given 
by (superscripts indicate which 'cat' the coordinates refer 
to) 



V 1 I = M 

Pt+i 



q p\+e{ql)+e'{ql,ql) 



and, 



1t+i 
Pt+i 



= M ( q l 
Pt 



e(q 2 ) + e'(qlq 2 ) 



(31) 



i r 2D 



Mil jl M[2 j2 Cjl j2 



(32) 



where 



C 



exp{[iK c N/(2n)} cos[(2ir / N^j 1 + j 2 )}} (33) 



and j 1 ,] 2 , 1 1 , 1 2 € {0, . . . , N — 1}. In the following calcu- 
lations we take K c = 1.0. 

We have studied the classical and quantum separability 
entropies. In all cases we have used the initial states 
given in terms of a Gaussian phase space distribution 
with dispersion equal to vh, and its quantum analogue, 
a coherent state on the torus. Both distributions are 
centered in the middle of the phase space, i.e., at (q,p) = 
(0.5,0.5). We have evolved 10 time steps (iterations) for 
all maps. 

In Fig. [T] we show the separability entropy [6| as a func- 
tion of time (in units of map steps) for the baker map 
(|25I26I) . We have used the same discretization number 
N = 2 9 for both, the classical and the quantum simu- 
lations and removed the effects of the torus periodicity 
on the Wigner distributions [l^. We have decomposed 
the 2-torus in coordinates q and p. While the classical 
separability entropy saturates to a value of order ~ 1 
bit, after an initial short transient growth, due to the 
exact solvability of its Liouville evolution [19| , and drops 
back towards zero for t ~ logiV due to coarse-graining 
of classical Liouvillean evolution, the Wigner separabil- 
ity entropy grows with time due to quantum interference 
and diffraction patterns produced by discontinuity of the 
map. The baker's map is therefore non generic, in that 
it exhibits very different classical and quantum results 
even within the Ehrenfest time range due to diffraction 
effects. 

In order to investigate a case with generic behavior 
we have taken the smooth (continuous) perturbed cat 
map (|29I30[) . In Fig. [2] we show the same quantities as 
in Fig. [T] Now there is quantum to classical correspon- 
dence, up to the Ehrenfest time t~E oc log TV. Note that, 
similarly to the case of baker map, the classical result also 
exhibits discretization time scale after which the classical 
distribution experiences a complexity reduction which is 
a numerical artifact due to finiteness of N. Indeed, due to 
coarse graining the Liouville distribution becomes homo- 
geneous in phase space, and this implies asymptotically 
vanishing separability entropy. Finally, and for compari- 
son purposes only, we show the quantum results without 
removing the effects of the ghost images that appear on 
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Time (map steps) 

FIG. 1: (Color online) Classical and quantum separability 
entropies behavior for the baker map. The (red) gray line 
with crosses corresponds to the classical result and the black 
line with squares to the quantum one. The quantum curve 
has been calculated using the Wigner function without ghost 
images 18!]. In both cases N = 2 9 . 



account of the torus periodicity (this is shown just for 
the lower value of N) . Wigner separability entropy with- 
out ghost images removed exhibits some deviations for 
shorter times (then, the fact that we have ghost peaks 
in Wigner distribution is more important than at later 
times), the removal of ghost images is therefore needed 
in order to obtain quantum to classical correspondence. 



Finally, we have tested the Wigner separability en- 
tropy and its relation to pure-state entanglement for 
two coupled perturbed cat maps (|31|32[) (in this case 
we have evaluated the Wigner distributions according to 
Ref. [13]). In practice we have decomposed the four- 
dimensional phase space in the coordinates (g 1 ,?? 1 ) and 
{<l 2 tP 2 ) corresponding to each of the two 'cats', respec- 
tively. We have also compared the results with the classi- 
cal separability entropy Q. It is worth mentioning that 
in order to obtain the classical complexity measure we 
have used the same phase space decomposition for the 
four dimensional Liouville distribution that we have pre- 
viously applied to the four-dimensional Wigner function. 
All these results are shown in Fig. [3] The agreement 
between the quantum measures is accurate within the 
machine precision, i.e. h[W] = 2E(\ip)) (in the figure 
we have rescaled von Neumann entropy by a factor of 
2), reflecting the match between singular value decom- 
positions for the density matrix and its Wigner distribu- 
tion. In this case we did not remove ghost images in the 
Wigner distributions since all the coherences measured 
by the separability entropy in these coordinate pairs are 
relevant to measure entanglement. Regarding the classi- 
cal separability entropy we find a reasonable agreement 
up to time ~ log TV where due to coarse graining h[p c ] 
starts to drop. It is important to underline that the clas- 
sical saturation is entirely due to the phase space coarse 
graining chosen for the numerical simulations, while the 
quantum one is an unavoidable phenomenon fixed by the 
finite size of h. Note that, as recently pointed out (2lj . 
the initial growth of entanglement can be reproduced in 
the semiclassical regime by purely classical computations. 
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Time (map steps) 

FIG. 2: (Color online) Same as in Fig. [TJ but now for the 
perturbed cat map. The black lines with squares correspond 
to quantum results and the gray (red) ones with crosses to 
classical ones. The light gray (green) dotted line corresponds 
to the quantum values using the Wigner distribution with 
ghost images which exhibits small deviations at shorter times. 
The lower (three) curves correspond to N — 2 11 while the 
upper (two) curves correspond to N = 2 13 . 




0123456789 10 
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FIG. 3: (Color online) Classical and quantum separability 
entropies behavior for two coupled cat maps (the gray (red) 
line with crosses and the black line with squares, respectively). 
Two times von Neumann entropy is also shown by means 
of a light gray (green) dotted line with triangles. We have 
used discretization number N = 2 6 in both the classical and 
quantum cases, for each map. 
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IV. CONCLUSIONS 

In this paper we have proposed a new measure of com- 
plexity of quantum states, the Wigner separability en- 
tropy. This quantity turns out to be equal to the oper- 
ator space entanglement entropy, it quantifies the mini- 
mal amount of computational resources required to sim- 
ulate the quantum dynamical evolution of a system by 
means of time-dependent density-matrix renormalization 
group Q. Moreover, due to its relation with the entropy 
of entanglement for pure states and to the existence of 
the analogous s-entropy of Liouville densities in classi- 
cal dynamics, the separability entropy emerges as an ex- 
tremely broad complexity quantifier in both the classical 
and quantum realms. 

With regard to a previously proposed phase-space 
quantum complexity indicator, that is, the number of 
Fourier harmonics of the Wigner function [H, 0|, the 
Wigner separability entropy has the advantage of being 
basis- independent. More importantly, numerical indica- 



tions [5] suggest that the Wigner separability entropy 
should be able to distinguish between quantum chaotic 
and quantum regular motion also for many-body sys- 
tems. Finally, we point out that the Wigner separability 
entropy is well defined also for mixed states and could 
therefore be used to quantify the complexity of decoher- 
ent quantum dynamics. 
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